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Dilute oil dispersions of fractal carbon black particles with attractive Van der Waals interac- 
tions display continuous shear thickening followed by shear thinning at high shear rates. The 
shear thickening transition occurs at 7c ~ 10^ — lO'^s"^ and is driven by hydrodynamic breakup 
of clusters. Pre-shearing dispersions at shear rates 7 > 7c produces enhanced-modulus gels where 
~ '^pre-shear directly proportioual to the residual stress in the gel measured at a fixed 

sample age. The observed data can be accounted for using a simple scaling model for the breakup 
of fractal clusters under shear stress. 

PACS numbers: 82.70.Dd,82.70.Gg,83.60.Rs 



Colloidal particles interacting with attractive poten- 
tials in a fluid form gels above some critical volume 
fraction, 0c, where (pc is a function of the interaction 
strength, U . The viscoelastic properties of the gels are 
determined by (/>, U and the topology of the network of 
paxticle contacts, characterized by a fractal dimension, 
df. Additionally, the properties are sensitive to the me- 
chanical history of the system. To obtain reproducible 
results. Theological studies of colloidal gels typically em- 
ploy a high rate pre-shear as a rejuvenation or initializa- 
tion step in which the mechanical history of the gel is 
effectively erased. The non-Newtonian flow curves asso- 
ciated with attractive colloidal suspensions typically dis- 
play monotonic shear thinning and as a result the mod- 
ulus established after rejuvenation is not overly sensi- 
tive to the exact pre-shear flow condition. Shear thicken- 
ing transitions are neither generally known nor predicted 
for systems where attractive interactions are sufficiently 
strong to induce flocculated gels [TJ |2] . Such transitions 
would be accompanied by drastic changes in the underly- 
ing fluid microstructure. As a result, the elasticity of gels 
would show a marked dependence on the nature of the 
pre-shear flow and its location on the flow curve relative 
to the shear thickening regime. 

Here we report on the observation of shear thickening 
for attractive colloidal particles in a simple Newtonian 
fluid. We flnd the elasticity of gels formed by pre-shearing 
above the shear thickening transition scales as a power 
law with the pre-shear stress, and is directly proportional 
to the internal stress in the sample, which we measure. 
We propose a scaling model which considers the depen- 
dence of the gel modulus on the cluster number density 
produced during pre-shear flow. The model shows good 
agreement with the experimental data. 

Measurements were made at 25° C on dispersions of 



' ^chinedum.osujl®yale.edu| Present Address: Department of Chem- 
ical Engineering, Yale University, New Haven CT 06511 



carbon black particles ranging from 2 to 8 wt.%. A 
non-polar, small-molecule fluid, tetradecane, was used 
as the suspending medium to avoid complicating influ- 
ences of electrostatics or adsorption of macromolecular 
solvent species onto the colloidal surface. The carbon 
black (Cabot Vulcan XC72R) exists as « 0.5/xm diame- 
ter fumed particles with fractal dimension — 2.2 and 
Pcarbon = 1-8 g/cm^ . Tetradecane {psoivent = 0.762 
g/cm^, rjsoivent = 2.8 mPa.s) was obtained from Aldrich 
Chem. Co.. Samples were prepared by vigorous mixing 
of the two components using a vortexer. Optical stud- 
ies were conducted using a Bohlin CS rheometer with a 
transparent base-plate through which samples were im- 
aged using a CCD camera. Rheological measurements 
were made in strain-control mode using a TA Instru- 
ments AR-G2 rheometer with cone-plate and double- wall 
Couette geometries. The fast hardware feedback of this 
instrument, when run in strain-control, enabled measure- 
ments of the internal stress within samples by applying 
continuously computed torques to maintain constant tool 
position. Samples were pre-sheared before each measure- 
ment to remove the effect of shear history using an em- 
pirically determined shear rate (100s~^) which provided 
consistent results [3] . Samples were checked using various 
geometries (standard and roughened) to ensure the ab- 
sence of rheological artifacts and phenomena such as wall 
slip and edge-fracture. They were also replicated using a 
pure strain-controlled instrument, ARES-LSI (Rheomet- 
rics). 

We measure the steady state viscosity as a function of 
shear rate using a cone-plate geometry. At low shear 
rates, between w lO*' and 10^s~^ optical studies us- 
ing the plate-plate geometry show a build-up of struc- 
ture in the sample via the formation of rolling, vorticity 
aligned cylindrical floes. Such shear-induced structures 
have been observed in other viscoelastic systems under 
flow such as thixotropic clay gels 4 , nanotube suspen- 
sions [5j and attractive emulsion droplets [6] . Viscometry 
of such "sticky" systems at low shear rates may be sus- 
ceptible to wall slip effects f7| and so caution is required 
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FIG. 1: Shear rate dependent viscosity for a series of weight 
fractions. + : 2%, > : 3%, Q : 4%, o : 5%, A : 6%, □ : 8% The 
dotted-Hne is the background viscosity of the solvent. Inset: 
Rate dependent shear stress. 

in interpreting any rheological data in this regime. At 
higher shear rates where homogeneous flows are readily 
achieved, 7 > 10s~^, Figure [T] there is modest shear 
thinning accompanied by the formation of densified clus- 
ters as shown in Figure |2^. Shear thinning persists until 
« 10^ — lO'^s"^ where there is a composition dependent 
continuous transition to a shear thickening flow. The 
flow curves exhibit a hysteresis loop on descending sweeps 
as is typical of thixotropic materials, and show good re- 
producibility on subsequent ascending-descending loops. 
Shear thickening has in fact been reported [5] , though not 
extensively considered, for a system consisting of carbon 
black particles in an adsorbing silicone oil. This was at- 
tributed to hydro-cluster formation due to the occurrence 
of stress overshooting at shear rates that produced shear 
thickening. In the present case, at high shear rates in 
the shear thickening regime, the system presents a finely 
dispersed microstructure in which the densified clusters 
produced at lower shear rates have been eroded. This 
results in a higher effective volume fraction of particles 
in the dispersion, and thus an enhanced viscosity. Figure 
[2]d. In contrast to the hard sphere case, shear thickening 
results not from hydro-cluster formation, but from the 
breakup of locally dense clusters of the fractal colloidal 
particles into less dense structures which increase viscous 
dissipation in the system. A negative thixotropy of ferric- 
oxide suspensions composed of acicular particles has also 
been attributed to a qualitatively similar mechanism [3]. 

Based on the dispersed particle size, we estimate the 
Peclet number, Pe — « 10^ — 10^, confirming that 
the high shear rate regime is dominated by hydrody- 
namic forces. Similarly, the particle Reynold's number 
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Re = ^^^^^ « 10~ — 10^"', indicating that viscous effects 
dominate and contributions from inertial forces should 
be minimal. Under similar conditions in Stokesian dy- 
namics simulations of hard spheres |101 lllj a transition 
is observed from positive to negative normal stresses on 
shear thickening due to hydrocluster formation at high 
Pe, with |A^i| ~ 7^ where Ni is the first normal stress. 
In the present system, inertia corrected Ni transitioned 
to slight net positive values at the onset of shear thick- 
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FIG. 2: Microstructure under shear in parallel plate geometry. 
Gap = lOOfim, = 3 wt.%, 1.5x1.5 mm area, (a) 7 — 133s~^ 
(b) 7 — 13303"^^. Illumination for (b) was several times higher 
than for (a) to obtain sufficient light transmission through the 
optically dense structure. Pixel binning was used to decrease 
the required exposure time. 

ening T^, suggesting, along with the optical data, that 
hydrocluster formation is not in fact the mechanism in 
operation here. The critical shear rate at the onset of 
shear thickening, 7c, scales roughly as (f>^ and the critical 
stress iTc as cj)'^ with /3 roughly between -1.5 and -2, and 
v fa 1. From measurements performed using high vis- 
cosity basestock oil {rj « 33 mPa.s), we observe that the 
critical shear rate is inversely proportional to the solvent 
viscosity, suggesting that shear thickening in this system 
is controlled by a critical stress, rather than by a critical 
rate. 

We use dynamic measurements to study the depen- 
dence of gel modulus on pre-shear flow and composition. 
Samples were sheared to remove the effects of flow his- 
tory and then pre-sheared at the rate of interest for 20 
minutes, more than sufficient to achieve a steady flow 
and constant viscosity. The system was then allowed 
to sit quiescently for 30 minutes. We measure the vis- 
coelastic storage and loss moduli, G' and G", in the 
linear regime, 7 = 0.05% at a fixed angular frequency 
u! = Irad/s., where 7 is the shear strain. The elastic 
modulus of gels formed by pre-shearing at a fixed rate 
in the shear thickening regime follows a typical power 
law G" 4'°', with exponent a — 3.5, as might be ex- 
pected for attractive systems p3_ . In the shear-thinning 
regime, before the thickening transition, the exponent is 
slightly higher, a = 4.2. The elasticity of thickened gels is 
strongly enhanced, typically by over one order of magni- 
tude, compared to gels produced by non shear-thickening 
flow. This suggests that in addition to volume fraction, 
(/), there is another parameter dependence involved in de- 
termining the elasticity of gels produced across a range 
of shear rates. The modulus follows the same qualitative 
dependence on shear rate as the viscosity, with a sud- 
den upturn at the critical rate or critical shear stress. A 
double-wall Couette cell was used to mitigate the effects 
of sedimentation at low shear stresses where samples did 
not form robust thickened gels, or were too dilute to be 
gravitationally stable for study in the cone-plate geome- 
try. Remarkably, as in Figure [sjinset), the dependence 
of the gel modulus on the pre-shear stress has a com- 
mon form across different compositions, (/», such that the 
data can be scaled onto a single curve using the criti- 
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cal stress, (Tc, and critical modulus, G^, as scaling pa- 
rameters. Data of shear thickened samples show a strik- 
ing dependence of the modulus on the pre-shear stress, 
G' ~ CT^'^~^, as shown in Figure [sjj main) . These samples 
which were studied in the cone-plate geometry at large 
pre-shear stresses where the resulting high modulus gels 
were not sensitive to sedimentation. There is a some- 
what weaker volume fraction dependence such that in 
the shear thickening regime, overall, the modulus of the 
colloidal gel is substantially influenced by the pre-shear 
stress applied to the system prior to gelation. 
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FIG. 3: Elasticity of shear thickened gels as a function of the 
pre-shear stress, using cone-plate geometry. ■ : 3%, A : 4%, 
• : 5%, ♦ : 6%, T : 8%. Inset: Data obtained spanning low 
and high stresses using a double- wall Couette cell. + : 2%, 
□ : 3%, O : 5% 

We interpret these results in terms of the dependence 
of the cluster size on shear stress. As the stress increases, 
densified clusters break, and assume a size that is set by 
the balance between the shear stress and the cohesive 
energy density of the cluster. 

The shear force on the cluster, Fgh^ar = 47r cri?^ where 
Rc is the radius of the fractal cluster. This force is spread 
across Ng particles in the plane of shear of the cluster 
where ~ _Rc' ^ where (i/ is the fractal dimension 
of the cluster. This force is balanced by the attractive 
interaction between particles, defined by the interaction 
strength U and the relevant length scale, a, the particle 
diameter. Thus we have 



which sets the dependence of the equilibrium cluster 
size on the shear stress. The elastic modulus scales di- 
rectly with the number density of clusters, i.e. G' ~ 
u U, the cohesive energy density, as in many disordered 
systems |14| . The number density, u of clusters of radius 

Rc scales as (p/R^^ so that overall, we expect the gel 
modulus to scale as 

G'-^CT^-"/ (2) 

Assuming a cluster fractal dimension of 1.8, as is found 
in diffusion limited aggregating systems then G' ~ a'^^'^4>. 



We plot this in Figure [4] and find good accord between 
our data and a slope of 1 as shown, over several decades. 
Our scaling is in rough agreement with empirical data 
on flocculated systems where the number of particles in a 
floe scales roughly as iVc ~ cr"^ Q^J. Although the simple 
assumptions of the scaling argument provide satisfactory 
results, the choice for df should be rationalized by more 
detailed measurements than we have so far made. 




FIG. 4: Elasticity of shear thickened gels as a function of 
the pre-shear stress re-scaled according to equation ([2| with 
df = 1.8. ■ : 3%, ▲ : 4%, •: 5%, ♦ : 6%, T : 8%. The solid 
line is a visual guide with a slope of 1. 

The abrupt cessation of shear in a freely flowing sys- 
tem above 0c results in a sudden entry into the gel state. 
The system loses mobility as particles experience sticky 
contacts with neighbors leading to the formation of struc- 
tures that rapidly span space, forming an elastic medium. 
This rapid sol-gel transition should be characterized by 
the development of internal or residual stresses, Ci, coin- 
cident with the mechanical quench - the network struc- 
ture, deformed with respect to an equilibrium relaxed 
configuration, would exhibit stresses which relax grad- 
ually in time. These stresses are believed to drive the 
peculiar dynamics observed in the aging of some glassy 
systems [HI [17] . Since they result from a deformation of 
the structure, the internal stresses are expected to be pro- 
portional to the elastic modulus [TH] , though no direct ex- 
perimental evidence as yet confirms this. Previous work 
has examined the very short time {t < 0.1s) dynamics of 
stresses in weakly fiocculating systems subjected to stress 
jumps |19j . Here, we followed ai{t) in our system from 
t K, Is after cessation of flow to t « 30 mins. and find a 
weak power law decay over this period. For shear thick- 
ened gels, the stress scales roughly as cr^ ~ t^"-^, inset 
Figure [Sj The internal stress acts counter to the direc- 
tion in which the sample was pre-sheared and is reversed 
on changing the pre-shear direction, indicating that the 
deformation of the clusters during flow (relative to the 
quiescent state) gives rise to the residual stress that is 
manifested by the rate quench and rapid gelation. Crit- 
ically, we find that the internal stress at a fixed sample 
age, i=30 mins., is directly proportional to the modulus 
of the system measured immediately thereafter, over sev- 
eral orders of magnitude of internal stress, produced by 
different pre-shear flows. Figure |5] The prototypical pic- 
ture of a slow diffusion limited gelation clearly does not 
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apply despite the relatively low volume fractions under 
study here. Gelation is unusually rapid and the system 
displays a well developed elastic modulus at the shortest 
times studied, ~ 1 second after cessation of shear. 




FIG. 5: Shear modulus scales linearly with residual stresses 
resulting from the quench into the gel state on cessation of 
flow. ■ : 3%, ▲ : 4%, •: 5%, ♦ : 6%, T : 8%. The solid hue is 
a visual guide with a slope of 1. Inset: Time dependence of 
internal stress, with open symbols for clarity. 

Shear thickening in hard sphere and repulsive Brown- 
ian systems has been rationalized using a force balance 
model that compares the shear force pushing particles to- 
gether to the diffusive flux that inhibits cluster formation 

Eol. 
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where a is the particle radius, h = r — 2a the distance 
between two particles, g(r) the particle distribution func- 
tion and fx the viscosity of the medium, replaceable by 
the suspension viscosity, ry, as a way of accommodating 
many body effects. In the system under study here, at 
low shear stresses, large non-Brownian clusters resulting 
from breakage of the gel at low shear rates have no dif- 
fusive barrier towards further clustering. They densify 
as shear stresses force them into contact, reducing their 



effective volume fraction in the system. Above a critical 
stress however, the cohesive energy density of the clus- 
ters is overcome and they start to disintegrate with a 
attendant rise in viscosity. 

The observation of shear thickening in this system 
of attractive particles was surprising. Indeed, previous 
careful examinations of flocculated systems have not re- 
sulted in observation of shear thickening in fractal j2] 
and "solid" particulate systems p]. A particularly sys- 
tematic examination of the role of interaction strength 
using depletion forces pLj shows that shear thickening di- 
minished and eventually disappeared on addition of sub- 
stantial attractive interactions. It is clear that in the 
current system, the fractal nature of the particles, irre- 
versibly fused agglomerates of smaller particles produced 
by flame pyrolysis, plays a key role. It drives the in- 
crease in the effective volume fraction on cluster break- 
up. Shear thickening is thus a product of shear induced 
reduction of cluster size and increase of the effective vol- 
ume fraction of the particles in the dispersion due to the 
fractal particle structure. Correspondingly, the modu- 
lus of shear thickened dispersions is enhanced due to the 
increase in cluster number density. The scaling of the 
elastic modulus with pre-shear stress is satisfactorily ac- 
counted for by consideration of the equilibrium cluster 
size as a function of this stress. The direct proportion- 
ality between the elastic modulus and the internal stress 
resulting from a system quench on cessation of shear has 
been confirmed by direct measurement. The dynamics of 
these internal stresses are a topic for continued studies. 
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